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A STUDY OF LIQUID BORIC OXIDE PARTICLE GROWTH RATES IN A GAS STREAM 
FROM A SIMULATED JET-ENGINE COMBUSTOR 
By Paul C. Setze 


SUMMARY 

It was experimentally determined that the liquid boric oxide particles 
leaving a jet-engine combustor that is burning a boron-containing fuel will 

have diameters of 1X10"® to 2.0X10“® centimeter. For this size range the 
particle drag coefficient and heat-transfer coefficient are essentially 
infinite. These results may be applied to any boron-containing fuel. 

A study of boric oxide deposition mechanisms is included and sugges- 
tions for decreasing deposition rates given. 


INTRODUCTION 

Boron hydrides are being tested as aircraft fuels because their heat- 
ing values are as much as 65 percent higher than those of hydrocarbons. 
They differ from conventional fuels in that, at the usual jet-engine op- 
erating temperatures, one of the combustion products, boric oxide (BgOj) , 

appears as a viscous liquid. The liquid boric oxide can produce perform- 
ance losses by collecting on critical engine surfaces and "spoiling" the 
flow of combustion gases. 

In order to calculate the performance of the boron fuels, it is nec- 
essary to know the growth rate of the liquid particles in the exhaust 

gases (ref. l) . When the particle diameter is less than 2.2X10“® centi- 
meter, 3 percent or more of the heat of vaporization is in the form of 
surface energy. When the particle consists of a single molecule, the 
entire heat of vaporization (1994 Btu/lb, ref. 2) appears as this sur- 
face energy and is not transferred to the gas stream. However, if the 

particle diameter exceeds 3X10 -4 centimeters, the entire heat of vapori- 
zation of the boric oxide that was previously held by the surface energy 
of the particles is absorbed by the working fluid, but the particles are 
so large that they cease to be in temperature and velocity equilibrium 
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with the gas stream. Reference 1 shows that for pentahorane fuel (BgHg) 

at a flight Mach number of 3.25 and a combustion temperature of 2200° R 
a maximum theoretical gain of 250 seconds in engine impulse is attain- 
able if the liquid boric oxide particles are in the range of infinite 
drag and infinite heat transfer (diameters less than 3X10"^ cm) . There- 
fore, it would be desirable for the particle diameters to fall between 

2.2x10"® and 3X10" 4 centimeter where less than 3 percent of the latent 
heat of vaporization is held by surface energy and the particles have 
essentially infinite drag coefficient and heat- transfer coefficient. 

Additional losses can occur when the liquid boric oxide deposits on 
engine surfaces. If the surface temperature is below the freezing point 
of boric oxide, large obstructions to flow may result. The manner in 
which boric oxide deposits and the techniques for reducing deposition 
rates are strongly dependent on particle size. 

Observations previous to this time (ref. 3) have indicated that the 
liquid boric oxide is in the form of very small spherical droplets that 
grow in volume as a function of time. These observations were purely 
qualitative, however, and no accurate estimate of the particle size or 
rate of growth was made. Along with these very scant experimental ob- 
servations, an attempt (ref. 3) was made to correlate the growth of liq- 
uid boric oxide particles with some of the parameters in a previous theo- 
retical work on the coagulation of colloids by Smoluchowski (ref. 4). 

In addition to this, a number of theoretical and experimental papers have 
been published on the coagulation of aerosols, but these systems do not 
simulate the particle growth process in a jet-engine combustor. 

Because of the status of this field of particle growth, an experi- 
mental and theoretical study of particle growth was initiated at the NACA 
Lewis laboratory. The two specific objectives of this program were: 

(l) establishment of particle size so that accurate quantitative values 
of performance parameters for boron fuels can be calculated; and (2) de- 
termination of size and motion of liquid boric oxide particles so that 
the mechanism of deposition can be better understood. 

The range of variables considered in this investigation was repre- 
sentative of turbojet-engine operating conditions: pressure, 1 to 1.5 

atmospheres; inlet temperature, 580° R; outlet temperature, 1260° to 
2260° R; and particle growth time, 2.5 to 30.0 milliseconds. In all cases 
the combustor outlet temperature was at least 1800° R lower than the 
saturation temperature of the liquid boric oxide. The boric oxide vapor- 
pressure data are taken from reference 5. 

The theoretical work took the form of a general particle growth 
equation based on statistical collisions between particles. This equa- 
tion was then applied to the two particle size ranges of interest: 
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(l) particle diameters less than 2.2x10"^ centimeter, where considerable 
energy is in the form of surface energy; and ( 2 ) particle diameters greater 

than 2.2X10 centimeter, where negligible energy is in the form of surface 
energy. 

The experimental program was performed in a 4- inch connected-pipe 
facility using pentaborane (B 5 H 9 ) fuel. At a predetermined station down- 
stream of the fuel injector, a portion of the exhaust gas entered a sam- 
pling probe, where it was rapidly quenched with liquid nitrogen to a tem- 
perature well below the freezing temperature of boric oxide. The now 
solid boric oxide particles were collected and analyzed for size, shape, 
and composition. 

To show the applicability of the results obtained in the 4- inch 
connected-pipe facility to turbojet engine applications, a single test 
was conducted using the single J47 combustor of reference 6 . 

The results of both the theoretical and experimental program are 
reported herein. 


SYMBOLS 

The following symbols are used in this report: 
A area, sq cm 

C concentration, g/cc 

Cp specific heat, cal/(g)(°K) 

D diameter, cm 

E energy, ergs 

f/a fuel-air ratio 

(f/a)* stoichiometric fuel-air ratio 

K diffusion coefficient, sq cm/sec 

k constant defined by equation ( 22 ) 

k' constant defined by equation (14) 

Z length, cm 

M molecular weight, g/mole 
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m mass of material diffused, g/sec 
N Avogadro's number, 6 . 02 X 10 ^ 

N c number of collisions 

n total number of particles 

n„, number of moles 


np number of particles per unit volume of gas 

n^ number of particles per unit mass of particles 


p pressure, dynes/sq cm 

R universal gas constant, dyne-cm/(mole) (°K) 

T absolute temperature, °K 
u velocity, cm/sec 


v volume, cc 


w mass, g 
X mole fraction 


x weight fraction 

6 diffusion boundary- layer thickness, cm 

X latent heat of vaporization, cal/g 
K molecular velocity, cm/sec 

p density, g/cc 

cr surface tension, dyne /cm 

t time, sec 

<P equivalence ratio, dimensionless 
Subscripts : 

F flame 



WACA RM E55I20a 


5 


f final 
g gaseous phase 
L - liquid phase 
0 at time zero 
p particle 
s mainstream 
T total 
w wall 
t at time t 

THEORY 

Theoretical Model and Preliminary Assumptions 

The system under consideration is that of an air stream flowing 
through a duct of constant cross-sectional area. A boron-containing fuel 
is added and the resultant fuel-air mixture burned. Gaseous boric oxide 
is formed in the flame zone. Additional air is added and cools the flow- 
ing stream to a temperatrue below the saturation temperature of boric 
oxide, (it is assumed that the boric oxide concentration is such that 
a negligible amount of boric oxide vapor will be present.) The boric 
oxide molecules can now collide and remain together forming a discrete 
liquid particle. These liquid particles will continue to collide and 
grow in size as a result of these collisions. 

The following assumptions regarding the system are made: 

(1) There is a range of conditions over which the collision of two 
or more liquid boric oxide particles will result in the formation of one 
larger particle (see refs. 7 and 8). 

(2) The particles are spherical. 

(3) The particles are evenly distributed across the duct cross 
section. 

(4=) At any one station in the combustor, the particle diameter range 
is sufficiently narrow so that an average particle diameter may be chosen 
that expresses the particle diameter term to be used in the growth equa- 
tion (see ref. 9). 
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(5) The particle may be treated as a large gas molecule. 

(6) The liquid droplets are treated as point particles. 


General Equations 


Assume that a particle of diameter Dp and velocity Up is mov- 
ing along a linear path. The moving particle now passes through a group 
of identical particles randomly distributed and at rest. The moving 
particle will collide with any other particle whose center lies within 
the distance Dp of the path along which it travels. During a unit of 

time x, the distance the particle travels is UpT. 

The number of collisions in time x is equal to the number of par- 

O 

tides whose centers lie within a cylinder of volume equal to itDpUpX. 

If n p is the number of particles per unit volume of gas, the number of 

* o 

collisions made by the moving particle in time x is equal to npitDpUpX. 


In reality all the particles are moving in random directions with 
average velocity Up. According to Loeb (ref. 10 ) , the velocity of any 

one particle relative to another is ^ Up. Therefore, the number of col- 
lisions per unit volume in unit time x is 


N C = 3 Itn p D P u P T 


( 1 ) 


and. the rate of collisions 


ai* c 4. 2 

dx = 3 


Assumption (l) states that the collision of two particles will re- 
sult in the formation of a single particle of volume equal to the sum 
of the volumes of the colliding particles. Therefore, the growth rate 
of a particle can be expressed as 


dx dx 


(3) 


Since a collision of two equal- size particles per unit time would double 
the particle volume, of three particles triple it, and so forth, 
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Therefore, 


dv-, a 2 

3^ - 1 ’VW'p 

The number of particles per unit volume is expressed as 


(5) 


V - 1 


v p v g 


Substituting equation (6) into equation (5) yields 


_ 4 rtJJ p u p v L 
dx = 3 v g 


( 6 ) 


(?) 


The findings of some investigators of Brownian motion are reviewed 
in reference 10. The results showed that the total average energy of 
particles undergoing Brownian motion was equal to the mean kinetic en- 
ergy of a gas molecule. In the application now under consideration, 
the particle diameter is less than the mean free path of the surround- 
ing fluid molecules and the particle motion will more closely approach 
pure molecular motion than Brownian motion. Therefore, the kinetic 
energy of the particles will be assumed to be equal to the kinetic en- 
ergy of a gas molecule, and the particle velocity equal to the velocity 
of an ideal gas molecule of corresponding mass. Reference 10 gives the 
velocity of an ideal gas molecule as (ref. 10) 



( 8 ) 


Wow 


and 


Therefore, 


M = WpN 


w p “ 6 ItD P p L 


„ 18RT 

? 3 

\ nD^o N y 


0 ) 


( 10 ) 


( 11 ) 
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Substituting equation (ll) into equation (7) gives 



For any set of inlet conditions 



Therefore, 



(Dp) 1 /* 


Since, 



then 

d p /2 A-?) 

Substituting equation (16) into equation (14) gives 



Setting limits for integration 



( 12 ) 


(l2a) 


(13) 


(14) 


(15) 


(16) 


(17) 


( 18 ) 
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which gives 


6 

5 




Substituting 



for 


v p 


yields 





Let 


(19) 


( 20 ) 


k 


•anr > 


= 1.591 k’ 


( 21 ) 


The general equation for the growth of a particle is 

D 5 / 2 - D 5 / 2 = kx 
p,x p,0 


( 22 ) 


where 


k = 



and T is in °K, and in grams per cubic centimeter. 


(24) 


The simple equation for the particle growth rate obtained is largely 
due to the initial assumptions. The validity of these assumptions will 
now be examined for the conditions treated in this investigation. These 
conditions as previously stated are 


Pressure, atm lto 1.5 

Temperature, °R 1260 to 2260 

Growth time, millisec (secXlO - ^) 2.5 to 30 

Supercooling, °R 1800 to 2200 


The assumptions subject to most question are numbers (l), the col- 
lision of two or more liquid particles will result in the formation of a 
single larger particle; and (4), a mean diameter may be assigned and 
will be typical of the particle diameters at any combustor station. 
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In regard to assumption (l), reference 7 states that 100-percent 
collision effectiveness may be expected. The experimental data of ref- 
erence 8 also verify assumption (l). The change in particle size due 
to evaporation caused by the surface energy released upon a collision 
is discussed in a following section. 

In regard to assumption (4) the particle distribution curves of 
reference 9 show that in the aging of smokes not more than a 7 -percent 
deviation from the average particle diameter can be expected in turbu- 
lent air for any given set of conditions if all particles have the same 
terminal diameter. This process is very similar to the one under con- 
sideration here. Also, in the nozzles of supersonic and hypersonic wind 
tunnels, observations of condensation of air and water vapor have been 
made and a very narrow range of particle diameters at any one station 
observed; however, there are no published data giving exact values of 
particle diameters. These observations along with the experimental data 
of this report (given in a later section) tend to validate this 
assumption. 


Particles Having Brownian Motion and Large Surface Energy 

When a liquid surface is in contact with a gas, a force known as 
surface tension cr is required to maintain the liquid molecules at the 
liquid surface. The total energy E required to maintain this force 
for a given surface area A is equal to aA. 

The surface energy of 1 gram of n£ liquid drops with surface ten- 
sion cr and uniform diameter Dp is 

E = oAn ' (24) 

P P 


Since 


and 






( 25 ) 
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This quantity E is large for particle diameters approaching the 
molecular diameter, and equals the latent heat of vaporization X (1108 
cal/g or 1994 Btu/lb for boric oxide at the normal boiling point) for 
single molecule particles. 

Figure 1 shows the percentage of the latent heat of vaporization 
contained in the surface energy of a given weight of liquid boric oxide 
particles of constant diameter (calculated from eq. (25)). 

When two particles a and b collide and form a new particle c, 
the surface energy released is 

(E a + Eh) - E c = AE (26) 

oAg^ + ctAb - aA c = AE (27) 

A^ and since v a + v^ = v c , D c = 2 1 / 3 D a . There- 

AE = (2 - 2 2 / 3 )A a a (28) 

The fraction of the initial surface energy released upon the collision 
of two particles of equal diameter is 


For example, if A a = 
fore. 


AE _ (2 - 2 2 / 3 ) 
2E a = 2 


0.206 


(29) 


The quantity of heat given up upon a collision of two particles as a 
function of particle diameter is shown in figure 2. This energy unless 
rejected to the surroundings will raise the temperature of the particle. 
If the temperature is raised above the saturation temperature, evapora- 
tion will reduce the size of the combined particle. (See CONCLUDING 
REMARKS.) 

Assume 1 pound of liquid boric oxide particles, each 5X10 centi- 
meter in diameter. If each particle collides with another particle, the 
heat released is calculated as follows: 

From figure 1, 18 percent of X is held by surface forces. There- 
fore, the heat release is 


0.18X1994X0.206 = 72 Btu 
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and the temperature rise 


(c p = 0.213 Btu/(°F)(lb) 

The rate at which this energy will be dissipated to the surround- 
ings will be proportional to the number of collisions the particle makes 
with gas molecules before it collides with another liquid particle and 
releases more energy. For a pentaborane equivalence ratio of 0.20 at 1 
atmosphere pressure and a particle diameter of 5X10“® centimeter, the 
time between liquid particle collisions from equation (22) is 1.62X10”H 
second. The number of collisions the original particle makes with gas 
molecules before colliding with another liquid particle is 2.9X102. If 
it is assumed that each colliding molecule comes to temperature equilib- 
rium with the liquid particle, this number of collisions is sufficient 
to reduce the temperature of the liquid particle enough so that negli- 
gible evaporation would occur. Calculations of the foregoing type indi- 
cate that where the degree of supercooling is at least 700° R it is un- 
necessary to consider surface energy at all and equation (22) may be used 
to express the particle growth rate in the region. If the supercooling 
is less than 700° R, it is possible that no condensation will occur; and 
in any event, a stepwise solution of equation (22) in conjunction with 
the result of equation (29) must be used to obtain particle growth rates. 

From this method of calculating stepwise particle growth, the time 
required for a particle to grow to a diameter of 2.2X10“® centimeter is 
less than 1X10“® second for 1800° R supercooling. The value of 2.2XL0"® 
centimeter is the minimum diameter where surface energy may be neglected 
as calculated by the method of reference 11. Growth times of the order 

of 1X10“6 second are of interest in the expansion of boron-fuel combus- 
tion products through a nozzle at high combustion temperatures (above 
3840° R unless the combustion pressure is greater than 1.5 atm, at which 
point temperatures as low as 3460° R should be considered) . Therefore, 
for the conditions under study in this report (maximum temperature, 

2260° R) the time required for growth to a diameter of 2.2X10”® centi- 
meter will be neglected. 

In reference 11 an analysis of the condensation of water vapor is 
presented. The time required for the liquid water particles to grow to 
a diameter where surface forces axe negligible (for water, 1X10“^ cm) is 
calculated to be less than 1X10"® second, which is similar to the value 
calculated for boric oxide. 


is 


72 

0.213 

rO 


= 339 u R 


or cal/( C)(g) from ref. 2). 
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The preceding analysis assumed spherical particles, which is not 
the case if only a small number of molecules make up a single particle. 
However, the effect of particles of this type is considered small be- 
cause of their short life and will be neglected. 


. Particles Having Negligible Surface Energy 

In the diameter range where surface energy may be neglected 
(Dp > 2. 2X10" 7 cm), equation (22) may be used to describe the particle 
growth directly. 


EXPERIMENTAL APPARATUS AND PROCEDURE 
General 

The experimental work was performed in a 4-inch connected-pipe 
facility supplied with metered combustion air from the laboratory central 
air system (fig. 3) . A perforated plate (fig. 4) blocked 70 percent of 
the pipe area l/4 inch upstream of the plane of fuel injection. A high 
pressure drop (Ap/p = 0.27 at an inlet-air velocity of 175 ft/sec) 
across the perforated plate served to achieve rapid mixing in the com- 
bustion zone and a flat temperature and velocity profile. Fuel was in- 
jected normal to the air flow from three points in each of six arms ex- 
tending through holes in the perforated plate (see figs. 4 and 5). 
Ignition was supplied by means of a spark plug 8 inches downstream of 
the perforated plate. 

The 18-inch-long section of pipe immediately downstream of the per- 
forated plate was replaced with a Lucite tube |4-in. I.D. by 4^- in . O.D.j 

for three preliminary tests. The purpose of these tests was to deter- 
mine temperature and velocity profiles and to observe visually and photo- 
graphically the pentaborane flame. Motion pictures of the flame were 
taken through the Lucite at a camera speed of 128 frames per second and 
a lens setting of f/ll. 

The fuel system (fig. 6) was designed for use with boron hydrides 
and provided for remote operation, flushing, and cleaning of fuel lines 
and injector. After each run, the injector was immediately flushed out 
with helium to prevent decomposition of the pentaborane fuel. The lines 
were then flushed with JP-4 fuel to remove excess pentaborane, which was 
followed by a 50 percent methanol - 50 percent acetone mixture tha t re- 
acted with any remaining pentaborane to form nontoxic noncombustible 
products. Finally the lines were again flushed with JP-4 fuel to remove 
all traces of fuel or fuel residue. Fuel flow was measured with a 
rotating- vane -type flowmeter and recorded on a self -balancing strip- 
chart recorder. 
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Temperature, static pressure, and total pressure were measured and 
recorded at three equidistant stations between the perforated plate and 
the sampling probe. These measurements were used to determine combus- 
tion efficiencies and temperature and velocity profiles. Combustion ef- 
ficiencies were computed by the method of reference 12. 


Particle Sampling Probe and Quench Technique 

A probe for collecting and freezing the liquid boric oxide combus- 
tion product was placed at one of three distances downstream of the per- 
forated plate: 36.5 inches, 53 inches, or 128.5 inches. The probe de- 

tail is shown in figure 7. Liquid nitrogen was supplied from a 100- liter 
Dewar flask under helium pressure and injected through the inner wall of 
the probe at a point 3 inches from the probe inlet. Gaseous nitrogen was 
also injected at this point and atomized the liquid nitrogen. The quanti- 
ty of liquid nitrogen injected was not measured, but was sufficient to 
drop the temperature of the collected stream to 250°±50° F. At this tem- 
perature the boric oxide in the exhaust stream is a solid and the other 
components gaseous. Temperature measurements showed that the quenching 
of the collected stream was virtually complete in a maximum distance of 
3 inches after injection of liquid nitrogen. 


Particle Collection System 

A 2-inch-diameter rubber hose carried the collected stream from the 
sampling probe to the particle collector. An electronic particle pre- 
cipitator (see fig. 3) was used to remove the solid particles from the 
gas stream. According to the manufacturer's specifications the collec- 
tion efficiency of the particle precipitator is 87 percent for particles 
1X10" ^ centimeter in diameter, and 99 percent for particles 9X10"^ cen- 
timeter in diameter. 

The collection medium was l/l6-inch Fiberglas matting, and was re- 
placed after each run. Figure 8 shows the appearance of the Fiberglas 
mat with its mounting frame after a typical run. The collected stream 
with solid particles removed was then returned to the main combustion 
gas stream. A butterfly valve located in the return line was used to 
control the stream velocity entering the probe. 


General Description of Topical Test 

Before each test the filter was loaded with new collection media, 
the fuel-system pressure checked and flushed with helium, and the liquid 
nitrogen lines hooked up. After air flow and inlet-air temperature were 
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established, the liquid-nitrogen flask was pressurized and 3 minutes 
allowed for the liquid nitrogen lines to. cool. The fuel bottle was 
pressurized and the fuel throttle set to give the desired fuel-air 
ratio. The fuel was turned on and once a flame was established (as 
evidenced by total pressure and temperature rise) 5 seconds were allowed 
for. instruments to come to equilibrium. The instrument recorders were 
then turned on and after 12 seconds (one instrument cycle) the fuel was 
turned off and the lines purged as described earlier. 


Particle Analysis 

The Fiberglas mat was removed from the particle precipitator and a 
small portion washed with n-heptane (C 5 H 12 ) to remove the boric oxide 
particles. Boric oxide is insoluble in n-heptane. A portion of the 
resulting suspension of boric oxide particles in n-heptane was deposited 
on electron microscope grids. The n-heptane was evaporated, leaving the 
boric oxide particles. A minimum of three electron-photomicrographs was 
taken of samples from each run. Examples of the electron-photomicrographs 
are shown in figure 9. The magnification of the electron microscope was 
calibrated and the size of each particle in the photographs measured. 
Approximately 50 particles were measured for each run. From these meas- 
urements an average particle diameter was determined by plotting ob- 
served diameter against percentage of total number of particles having 
that diameter (see fig. 10) . The particle agglomerates shown were formed 
in the particle precipitator and the heptane suspension. Careful visual 
observation of the agglomerates under the electron microscope showed they 
were made up of a number of smaller particles of approximately spherical 
shape, some of which are outlined in figure 11. Fifty measurements were 
considered sufficient for the analysis because of the very narrow range 
of sizes encountered. The probable error in average particle diameter 
was calculated to be ±11.4 percent for run 6, which gave the largest 
spread of measured particle diameters. The error in actual measurement 
of the particle diameters was within ±3 percent. 


Calculation of Liquid Particle Residence Time 

The flame length was determined by temperature measurements and 
photography. It was assumed that the pentaborane enters the combustor 
as a vaporizing liquid droplet surrounded by air and burns near the 
stoichiometric flame temperature. The. boric oxide thus formed is in 
the vapor state and is later condensed by the rapid mixing with addi- 
tional air. The point where the condensation of the boric oxide occurs 


I 
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was assumed to coincide with the point where the flame becomes non- 
luminous. With these assumptions the time available for particle growth 
t can be expressed by 



(30) 


where Z T is the distance from the perforated plate to the liquid nitro- 
gen quench station; Zp is the photographically determined flame length; 
and u s is the mean stream velocity. In some tests, the photographical- 
ly determined flame length was compared with the values indicated by 
thermocouple measurements. 


Range of Experimental Program 

The experimental program covered a pentaborane equi valence -ratio 
(actual fuel-air ratio divided by stoichiometric fuel-air ratio) range 
of 0.074 to 0.22 and a combustor static-pressure range of 1 to 1.4 at- 
mospheres. The inlet-air temperature was held constant at 580° R and 
the inlet-air velocity varied from 130 to 250 feet per second. The 
total time available for particle growth varied from 2.5 to 30 milli- 
seconds. 


RESULTS 

Comparison of Flame Lengths Determined from 
Photographs and from Thermocouples 

From the photographic studies made with the replaceable Lucite tube, 
a flame length of 16 to 18 inches was observed over an equivalence-ratio 
range of 0.10 to 0.24. The use of visible flame length as a criterion 
for determining the length of the combustion zone was partially con- 
firmed by making thermocouple measurements during representative runs . 

For example, in preliminary test 2 the temperature was 1800° F 22 inches 
downstream of the fuel injector while a thermocouple 16 inches down- 
stream of the fuel injector was burned off, which indicates that the 
combustion reaction had not been completed at that point. A flame length 
of 18 inches was observed from motion pictures. Therefore, when using 
equation (35) for determining time available for particle growth, a value 
of 18 inches was used for Z^. A photograph of the pentaborane flame in 

the Lucite section is shown in figure 11. 
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Combustion Efficiency 

A combustion efficiency of 100 percent was achieved during each of 
the tests. Also a chemical analysis of the collected samples showed 
them to be 100 percent boric oxide. 


Temperature and Velocity Profile 

Figure 12 shows the temperature and velocity profiles for run 2 , 
which was typical for each of the runs. For any one run the maximum 
temperature -profile spread from the mean temperature was 5 percent or 
less, and the maximum velocity profile spread from the mean velocity 10 
percent or less at a station 22 inches downstream of the perforated 
plate. The highest velocity occurred about l/2 inch from the center 
and was due to the penetration of the jets through the holes in the 
perforated plate. These jets decayed very rapidly from this point on, 
and were not noticeable at a point 30 inches downstream of the fuel 
injector. 


Particle Size and Growth Rate 

A summary of the experimental data appears in table I. 

Over the entire range of experimental conditions covered, the meas- 
ured average particle diameters fell between 1 . 10X10“ ^ and 1.46X10”^ 
centimeter, and the particles were very nearly spherical in shape. The 
experimentally measured particle diameters are shown plotted against 
particle diameter calculated from equation (22) for the same conditions 
and residence time in figure 13. The average deviation of the measured 
diameter from the calculated diameter was 27 percent, which gives a 
probable error of ,±34.1 percent. If the sign of the deviation is taken 
into account (greater or less than theoretical), the arithmetic average 
deviation is 16.4 percent. 


Particle Size in Turbojet Combustor 

The perforated-plate -type combustor and connected pipe used in the 
experimental work very closely approximate a ram-jet-type combustor; 
therefore, the experimental results are very easily applied to ram- jet 
calculations. To verify the experimental results for use on turbojet 
applications, a single run was made with the sampling probe at the out- 
let of the J47 tubular combustor reported in reference 6. Pentaborane 
fuel was burned at an equivalence ratio of 0.176. The reference veloc- 
ity was 103 feet per second and the inlet-air temperature 580° R. The 
boron oxide particle size was 1.23X10" 5 centimeter. This diameter 
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corresponds to a residence time of about 7 milliseconds. The exact 
residence time in a tubular combustor is unknown, but 7 milliseconds is 
a reasonable value. 


DISCUSSION AND APPLICATION OF RESULTS 
General 

It appears from both the theoretical and experimental data that the 
liquid boric oxide particles found in the combustion products of a jet 
engine burning a boron-containing fuel grow very rapidly to a diameter 
of about 1.0X10"^ centimeter. In the normal jet-engine residence time, 
the particles would never exceed a diameter of about 2.0X10-5 centimeter 
unless allowed to contact a surface where larger particles would be torn 
from the resultant film. This is illustrated in figures 14 and 15. 

The experimental data do not constitute an extremely critical test 
of the theory. The rather narrow range of particle diameters encountered 
probably accounts for the good data correlation obtained. Difficulties 
encountered in making accurate particle- size measurements probably ac- 
count for some of the deviation of the experimental results from the 
theory. However, equation (22) will give agreement within ±45 percent 
with the experimental results when calculating liquid boric oxide par- 
ticle diameters. 

The very narrow range of particle sizes encountered might indicate 
that condensation occurred in the collection probe. However, it was 
felt that this was not the case for the following reasons: 

(1) 100 Percent combustion efficiency was measured from the temper- 
ature rise across the combustor. Because of the high heat of vaporiza- 
tion of boric oxide, the temperature rise would have been about 18 per- 
cent less if the boric oxide had not been in the liquid state. 

(2) The degree of supercooling was sufficient for rapid condensation 
to occur immediately downstream of the flame zone. 

(3) The very short time available for particle growth allowed after 
injection of the liquid nitrogen (before fusion occurred) would make it 
very unlikely that the particles would grow to the diameters encountered. 

It is fortunate, from the standpoint of over-all engine efficiency, 
that the particle diameters are in the size range of IX 10“ 5 to 1.5X10-5 
centimeter. These particle diameters are large enough so that the sur- 
face energy is only about 1 percent or less of the latent heat of vapor- 
ization and yet small enough that the drag coefficients and heat-transfer 
coefficients based on the particle diameter are essentially infinite. 
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Performance Calculations 

By using equation (22), a boric oxide particle diameter may be cal- 
culated for conditions in the range of interest. This equation is ex- 
pressed in general terms, and may be used regardless of the fuel used 
(boron-hydrogen or boron-carbon-hydrogen) . This will make analytical 
evaluations of drag forces acting on the particle and heat-transfer 
rates to and from the particles (ref. l) possible. These calculations, 
in turn, will make more accurate analytical determinations of performance 
parameters attainable for all boron-containing fuels. 


Oxide Deposition on Engine Surfaces 

From the experimentally determined particle diameters it appears 
likely that most of the boric oxide found on engine surfaces is depos- 
ited by particle diffusion, or at least may be treated as a diffusional 
process. The rate of diffusion of a substance to a surface may be ex- 
pressed, in general, by Fick's law (ref. 13) 

S . K f (31) 

where K is a function of the type of diffusion occurring (turbulent or 
molecular) . 

Integration of equation (3l) from 0 to Z and from C s to C w gives 


m = K 



(32) 


Assume that all particles coming in contact with a wall remain on 
the wall, as would be the case with viscous liquid droplets. Therefore, 
the concentration of particles in the gas adjacent to the wall will ap- 
proach zero. Equation (32) becomes 

C„ - 0 (LA 

m = K A = K ~ (33) 

The particle concentration in the main gas stream C s is finite 
and essentially constant across any flow-channel cross section. A steep 
concentration gradient must then exist across a thin gas film adjacent 
to the wall, which is the diffusion boundary layer. 

The flow pattern of gases in a jet-engine combustor is extremely 
complex and not well enough defined to permit a calculation of the dif- 
fusion boundary-layer thickness 5. Also, the diffusion mechanism is 
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not defined, and, consequently, K is unknown. The simplest method of 
accurately determining 5 and K would be experimentally on carefully 
controlled surfaces in a combustor. Once 5 and K are known, deposi- 
tion rates on engine surfaces may be calculated from equation (33) by 
substituting 5 for l and the correct value of K. 

In view of the small particle sizes, it is now possible to suggest 
some methods that might eliminate or at least reduce engine deposits. 

It is seen from equation (33) that the rate of deposition on a sur- 
face will be inversely proportional to 5. Increasing S several times 
would decrease the deposition rate by a like amount. This increase in dif- 
fusion boundary-layer thickness may be accomplished by several methods: 

(l) introducing an air film along the surface on which deposition is oc- 
curring; (2) cooling this air film to freeze the particles before they hit 
the wall so that they might be easily swept off; or (3) introducing a 
mechanical device that would induce boundary- layer growth. 

The effect of a cold-air film on deposition rates is illustrated when 
the conventional combustor-liner deposition rates of references 14 and 15 
are compared with the deposition rate on the experimental wire-cloth com- 
bustor liner of reference 16. This comparison is shown in figure 16. 

The wire cloth allowed a continuous flow of secondary combustion air which 
filmed the liner surface, reducing boric oxide deposits 100 times. 

For particles formed by the growth process discussed herein, there is 
little possibility of any particle in the gas stream being large enough to 
contact a surface because of its own inertia. This could conceivably oc- 
cur, however, if particles of larger size were formed by some other means. 
One probable source of larger particles is as follows: With improper 

fuel-injector design some fuel is sprayed directly onto the walls of the 
combustor. This results in films of molten boric oxide (which might also 
be deposited by diffusion) on the walls. As these films flow past sharp 
corners or air-entry ports, large particles are torn from the film. To 
eliminate particles of this type, care must be taken in the design of en- 
gine components to eliminate surfaces that might generate films of molten 
boric oxide. 


CONCLUDING REMARKS 

The calculations which resulted in figures 1 and 2 were based on the 
assumptions that the surface tension of the infinite liquid is valid when 
applied to small liquid particles made up of a very few molecules, and 
that these particles are spherical in shape. Both of these assumptions 
are difficult to justify, but were made for lack of anything better. Con- 
sequently, the error involved in using figures 1 and 2 is probably large. 
However, these figures indicate a trend and are probably within an order 
of magnitude of the actual values. 

Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, September 26, 1955 
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APPENDIX - METHOD OF CALCULATING v B /v g FOR ANY 

FUEL CONTAINING BORON, CARBON, AND HYDROGEN 
AT COMBUSTION TEMPERATURES BELOW 2500° F 
Calculation of Vg/v g 


With 1 pound of inlet air as a basis, the total weight of oxygen # 
used by a combustion reaction at the stoichiometric fuel-air ratio (f/a) 
is 0.2329 pound. An oxygen balance for the combustion of a fuel con- 
taining boron, carbon, and hydrogen gives the following equation: 


0.2329 = 


x B (f'/ a )* 

3 

v '-‘ y 'ko nn 

x g(f/a) * 

on 4- 

x H 2 ( f / a )* 

m b J 

A . aOu • UL/ T 
4 

Me J 

A JL • \J\J T 

m h 2 _ 


|XtjX32.00 


(Al) 


Substituting the correct values for Mg, Mg, and Mg 2 (the molec- 
ular weights of boron, carbon, and hydrogen, respectively) and solving 
for (f/a)* gives 



1 

" 9.552 x B + 11.44 x c + 34.07 x R 


(A2) 


which gives the stoichiometric fuel-air ratio as a function of the 
weight fraction of boron, carbon, and hydrogen in the fuel (x £ , Xg, and 

X H 2 > • 


Equivalence ratio <P is defined as 

(f/a) 

actual 


cp = 


(f/a)* 


(A3) 


The total molar flow of exhaust gases per pound of inlet air n m ^ T 

(assuming all the boron oxide is in the liquid state and will be neg- 
lected) may be expressed by 


n^ = 0.02737 + 0.02272 


(fVc> ■+ 9 - 055 ^ 


(x H ) + (l - <P) (0.007278) 
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If the combustion temperature T c and combustion static pressure 
p c are known, the volumetric flow of the exhaust gases Vg per pound 
of inlet air can be calculated, using 



g 


°° 750 Pm.lffc 
Pc 


(A5) 


The volume of liquid in the exhaust gases vl per pound of inlet 
air Isx calculated from 



(A 6 ) 


The value of (density of boric oxide, glass) is taken from figure 17. 

.. From 'equations (A5) and (A 6 ), the volume concentration of boron ox- 
ide v i/ v g is easily calculated for any boron-carbon-hydrogen fuel. 
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TABLE I. - EXPERIMENTAL DATA BASED ON 1 GRAM OF INLET AIR PER SECOND 


Fuel Inlet-air Equiv- Combus- Combustioi 

flow, reference alence tion pressure, 

g/sec velocity, ratio, temper- p, 

ft/sec <P ature, atm 

T, 


-Combustion Weight Volume 


|0. 01305 
.01561 
.008389 
.007938 
.005661 
.008198 
.01065 
.01075 
.01040 
.01233 
.009947 
.01683 



0.1708 

.2043 

.1098 

.1039 

.0741 

.1073 

.1349 

.1407 

.1361 

.1614 

.1302 

.2203 



1.058 

1.233 

1.303 

1.419 

1.416 

1.169 

1.188 

1.128 

1.138 


flow 
of boric 
oxide , 
g/sec 


3.602X10' 

4.308 

2.315 

2.191 

1.562 

2.263 

2.939 

2.967 

2.870 

3.403 

2.745 

4.645 


flow 
of boric 
oxide, 

V L> 

Q.Q.J sec 

2.345X10' 

2.837 

1.476 

1.393 

.9842 

1.441 

1.892 

1.911 

1.846 

2.208 

1.761 

3.076 


Run Particle Total 

diameter growth 

calculated time, 

from equa- t , 

tion (22) 


1.025X10' 

1.005 

.842 

1.620 

1.480 

1.223 

1.800 

1.741 

1.528 

1.522 

1.458 

.790 


5.384 

5.763 

30.12 

24.57 

21.09 

22.71 

20.54 

19.06 

16.00 

17.94 

2.493 


Final 

Combus- 

Boric ox- 

Experi- 

particle 

tor 

ide volume 

mental 

diameter, 

exit 

concentra- 

d p 

D P,t> 

veloc- 

tion , 

Calcu- 

cm 

ity, 

v T /v , 
L' g’ 

lated 


u s> 

cc/cc 

D p 


ft/sec 



1.17X10 -5 

458 

9.650X10" 6 

1.141 

1.37 

541 

10.23 

1.363 

1.20 

506 

7.704 

1.425 

1.32 

306 

7.567 

.815 

1.46 

375 

7.130 

.986 

1.42 

436 

9.518 

1.161 

1.15 

405 

11.91 

.639 

1.10 

448 

11.93 

.632 

1.25 

483 

9.516 

.821 

1.35 

576 

10.63 

.887 

1.16 

513 

9.054 

.796 

1.37 

619 

11.75 

1.730 


perforated 
plate to 
quench 
station, 
in. 
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Figure 2. - Percent of latent heat of vaporization evolved upon collision 
of two identical particles. Temperature, 1980° R. 
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Figure 4. - Perforated plate and fuel injector after 70 seconds of operation (preliminary 
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Figure 5. - Fuel injector. 
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A Liquid nitrogen for injection 
B Liquid nitrogen for cooling 
C Gaseous nitrogen for atomization 
D Inner wall static -pres sure tap 
E Outer wall static -pres sure tap 



Liquid nitrogen 
^3 Cooling nitrogen (liquid) 
Atomizing nitrogen (gaseous) 


Figure 7* - Detail of sampling probe. 
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Figure 8. - Fiberglas collection media after typical run (run 9). 
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Figure 9. - Typical boric oxide particles 
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Figure 10. - Typical particle size distributions. 


MCA EM E55I20a 


35 



C-39964 



1 " 

2 


9" 




Figure 11. - Pentaborane flame through Lucite test section. Inlet- 
air velocity, 105 feet per second; equivalence ratio, 0.22. 



Temperature, °R Velocity, ft/sec 
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Figure 12. - Temperature and velocity profiles 22 
inches from perforated plate. Run 2; equivalence 
ratio, 0.20. 




Theoretical particle diameter after time T as calculated 

from eq. (22), cm 


MCA BM E55I20a 


37 



Experimentally measured particle diameter after time cm 

Figure 13. - Comparison of experimental results with 
equation (22) . 
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Figure 14. - Theoretical particle growth rates 




Deposition rate, g/sec 
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q.j no/qi ^q-fsusp (ss-spS) appxo opaog 


NACA - Langley Field, Va. 


Figure 17. - Density of boric oxide as a function of temperature 



